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For the theoretical prediction of cross-section fluctuations in chaotic scattering, the cross-section autocorrelation 
function is needed. That function is not known analytically. Using experimental data and numerical simulations, we 
show that an analytical approximation to the cross-section autocorrelation function can be obtained with the help of 
expressions first derived by Davis and Boose. Given the values of the average S-matrix elements and the mean level 
density of the scattering system, one can then reliably predict cross-section fluctuations. 



1. Purpose 

Quantum chaotic scattering is an ubiquitous phenomenon. 
It occurs, for instance, in nuclear physics in electron 
transport through disordered mesoscopic samples @], and 
in microwave billiards Q. In all cases, the cross section 
displays random fluctuations versus energy or frequency. 
These are due to the random features of the underlying 
resonances. With d the average resonance spacing and T 
the average width, data on cross-section fluctuations exist 
for the entire range of the parameter T/d, from the regime 
of isolated resonances (r <C d) to that of strongly overlap- 
ping resonances (r » d) . The analysis of the data focuses 
on the value of the average cross section and on quantities 
that characterize the cross-section fluctuations. These are 
the variance of the cross section and certain correlation 
functions. For the analysis, one needs theoretical expres- 
sions for these quantities. These should be generic and 
only use a minimum of adjustable parameters. 

The generic theoretical treatment of chaotic scattering 
employs a combination of scattering theory and random- 
matrix theory [l| and uses as input the values of d and of 
the energy-averaged elements S of the scattering matrix S. 
Analytical results exist for the 5-matrix autocorrelation 
and cross-correlation functions [l[ (including the value of 
the average cross section) and for the third and fourth 
moments of the 5-matrix Because of the complexity 

of the problem, analytical results for higher moments of 
the cross section or for cross-section correlation functions 
cannot be expected in the foreseeable future. 

The present paper aims at filling that gap. We com- 
bine the available analytical information [J [j, [f| , results 
of computer simulations, and of experimental work on mi- 
crowave billiards @, 0] to study the cross-section autocor- 
relation function for all values of T/d. In particular, we 



address the following questions, (i) For which values of 
T/d and with which accuracy can the cross-section auto- 
correlation function be predicted in terms of the S-matrix 
autocorrelation function? (ii) Which analytical alterna- 
tives exist should that approach fail? 

2. Framework 

We consider chaotic scattering in a time-reversal invari- 
ant system described by a unitary and symmetric scat- 
tering matrix S a b{E). Here a, b — 1,...,A denote the 
channels and E the energy (or, in the case of microwave 
billiards, the frequency). The number A of channels may 
range from unity to a large number, A> 1. Chaotic scat- 
tering is modeled by writing the S'-matrix in the form 7] 



Sab(E) 



W aa D~l{E)W hv 



where 
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The real and symmetric Hamiltonian matrix H has di- 
mension N ^> 1 and describes the dynamics of N res- 
onances labeled by Greek letters. These are coupled to 
the channels by the real matrix elements W afJ _. Chaos is 
taken into account by choosing H as a member of the 
Gaussian orthogonal ensemble (GOE) of random matri- 
ces Thus, the elements of H are Gaussian random 
variables with zero mean values and second moments given 

8„ 



by H av H pcr 



= {\ 2 /N)[5 ap u ;n - 



Sua + S pcr 5 V p\. Here and in the 
sequel, the overbar denotes the average over the ensemble. 
The parameter A determines (or is determined by) the av- 
erage level spacing d of the N resonances. It is convenient 
to decompose S(E) into an average and a fluctuating part, 
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Sab(E) — S aa S a b + Sab(E) 



(3) 
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where the average S-matrix S is assumed to be diagonal. 
The values of the diagonal elements S aa serve as input pa- 
rameters for the statistical model and are assumed to be 
known. That is the typical case: In nuclei, S aa is given 
in terms of the optical model of elastic scattering, in mi- 
crowave billiards S aa is determined by the running average 
over a measured spectrum 3, 6]. In rare cases, the average 
S'-matrix may not be diagonal. By an orthogonal trans- 
formation in channel space, S can be reduced to diagonal 
form, see Refs. [5|, [8|. For simplicity we do not address 
that case. By the same transformation, the phases of the 
S'-matrix usually appearing as factors on the right-hand 
side of Eq. ([3]), can be removed. Both for the S-matrix 
model Eq. |T]) considered in the present work and the ex- 
perimental data the average S-matrix is real and diagonal. 
Starting from Eq. (JlJ, the S-matrix autocorrelation func- 
tion (or "two-point function" ) 



(4) 



has been calculated analytically [l| for TV 3> 1 and fixed A. 
The resulting expression depends only on the difference e 
of the two energy arguments, on the average level spacing 
d of the system, and on the transmission coefficients T a of 
all channels a defined by 



\S a a I 



(5) 



The transmission coefficients obey < T a < 1. These 
coefficients measure the unitarity deficit of the average S- 
matrix and give the probability with which the resonances 
take part in the reaction. This is seen by using the de- 
composition Eq. ([3]) and the definition Eq. §5§ to write the 
unitarity condition for S in the form 



E 

b 



(6) 



For T a = or |S aa | = 1, we have S^ b (E) = for all b, and 
the resonances are not reached from channel a. Conversely, 
J2b l^ab^)! 2 i s maximal for T a = 1 or S aa = (complete 
absorption of the incident flux in channel a by resonance 
formation). The transmission coefficients T a determine the 
average width T of the resonances. An approximation for 
F is the "Weisskopf estimate" 



r 



|5> 



(7) 



The case of stron gly overlapping resonances r 3> d ("Er- 
icson regime" foL 1 1 Ql. 11 1) occurs for ^2 a T a 3> 1: The 
number A of channels must be large and most of the in- 
dividual transmission coefficients T a must not be small. 
Conversely, the case of nearly isolated resonances T <C d is 
realized when A is of order unity or when A is large but all 
T a are small. The theory developed in Ref. [l( and used 
in Refs. 0, [B| applies to all values of T/d. Equation ([7]) is 
exact in the Ericson regime and fairly reliable elsewhere. 



Under omission of kinematical factors the cross section 
in nuclear physics, the conductance in electron transport 
and the transmitted power in microwave billiards are all 
given by |S ab (-E)| 2 or by a sum of such terms. For brevity 
we refer to \S ab (E)\ 2 as to the cross section. The average 
cross section |S ab (£)| 2 = \Sa^{E)\ 2 6 ab + \S^ b (E)\ 2 is given 
in terms of S aa and of C„? (0) and is, thus, known. Fluc- 
tuations of the cross section are measured in terms of the 
cross-section autocorrelation function 



C ab (e) = \S ab (E + e/2)\ 2 \S ab (E - e/2)| 2 - \S ab \ 



(8) 



That function is the object of central interest in the present 
paper. With the help of the decomposition Eq. §5§ we 
write 



C ab ( e ) = 2(5 Qb ^e<^S Qa 2 C( 2 )(£) 



S aa S*Z(E + e/2)\Sfi a (E-e/2)\* 
S^S^(E-e/2)\S^E + e/2)\' 

\S? ib (E + e/2W\S« ab (E-e/2)\i- 



\s R J 2 



(9) 



We have used that in the experiments and in the consid- 
ered S-matrix model (Eq. ([I])) S aa is real, that by defini- 
tion S% b {E) = and that S» b (^i)S« fc (£; 2 ) = for all a, b 
and all E%,E2- The last relation holds because all poles 
of S lie in the lower half of the complex energy plane. To 
determine C ab (s) we need to know the four-point function 



Cat\s) = \S%(E + e/2)\z\S* b (E-e/2W 



\s a ab \ 2 



(10) 



and, in the elastic case a = 6, also the three-point function 
Cf b \e) = S^(E + e/2)\S* b (E-e/2)\ 2 . (11) 

These functions are known analytically only for e = and 
in the Ericson regime (r>d), see below. 

To determine magnitude and e-dependence of C ab (e), 
we combine analytical results with numerical and exper- 
imental evidence as follows, (i) Analytical results: In 
Refs. 0, HI analytic expressions are given for two func- 
tions F^t ' (e) and F^ b (e) that look similar to but actually 

differ from Cj$(e) and (e), respectively. These are 
defined by 



F%(e) = [S*J(E + e/2)Y[S* ab (E-e/2 
F^(e) = S*J(E + e/2) [S* ab (E - s/2)} 2 



(12) 



We note that in C ab (e) and in C ab (e) the elements Sf b 
and S^ b carry pairwise the same energy arguments. This 



is not the case for F^ b \e) and F^'(e). However, F^ b '(0) 



(3), 



p(3). 



coincides with C$(0), and F^'{0) differs from C^(0) 

— a — 2 

only by the known term |S" b | 2 . We are going to show 
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that for e ^ it is possible to approximate C ab (e) in 



terms of F^ b {e), and under certain conditions C^ ; (e) in 

terms of F^ b (e). For the convenience of the reader we, 
therefore, give in the Appendix analytic expressions for 



-«(3), 



ee C$( e ). We 



F^\e) for n = 2,3,4, where F%>(e) 
briefly show how the threefold integrals can be evaluated 
numerically to avoid the apparent singularities of the in- 
tegrand, (ii) Numerical results: For the numerical simu- 
lations we use Eqs. (UJ and ^ and fixed values for A, for 
the transmission coefficients T a , and for N as initial val- 
ues. Calculations were typically done for several 100 re- 
alizations to minimize statistical errors. The results agree 
very well with the available analytical results but go be- 
yond them, (iii) Data: The data stem from measurements 
of transmission and reflection amplitudes of microwaves 
in a flat cylindrical resonator made of copper and cou- 
pled to two antennas, see Refs. @, 0|. Microwave power 
was coupled into the resonator with the help of a vector 
network analyzer. The range of the excitation frequency 
was chosen such that only one vertical electric field mode 
is excited. Then the microwave cavity simulates a two- 
dimensional quantum billiard 12, Tj^. The resonator had 
the shape of a tilted stadium billiard whose classical dy- 
namics is chaotic. Transmission and reflection amplitudes 
correspond to complex S-matrix elements that are theo- 
retically modeled by Eqs. © and © @,|3|. 

3. Inelastic case (a y£ b) 

The inelastic case is simpler than the elastic one be- 
cause it involves only the function cjj^ (s) , see Eq. © . We 
begin with the Ericson regime T ^> d, see Refs. [9. I lot lll|. 
In Ref. [ll[ it was suggested and in Ref. [l5| it was shown 
that for T > cl the fluctuating S'-matrix elements S^ b 
possess a bivariate Gaussian distribution centered at zero. 
That fact implies that all higher moments and correlation 

(21 

functions can be computed from C ab (e) by way of Wick 
contraction. Hence C ah (e) = and 



C%(e) = \C%(e)\ 



(2), 



(13) 



other cases result from numerical simulations with, respec- 
tively, A = 32 and A = 52 equal transmission coefficients 
and imply, that Eq. (fTB)) holds for T > 3d. The ratio 

Ci2''(0)/I^'i2' ) (0)| 2 increases dramatically with decreasing 
T/d so that Eq. (fTB)) cannot be used much below r ss d. 



In the Ericson regime, the two-point function has the value [9|, 



C%(e) = (l + 5 ab ) 



T a T b 



(14) 



J2 c T c + 2iire/d ' 

Thus, for a ^ b and T ^> d, the cross-section autocorrela- 
tion function is known analytically. It has the shape of a 
Lorentzian with width T as given by Eq. (J7|). 

How far can we use Eq. (fTB")) outside the Ericson regime, 
i.e., for smaller values of T/d? Figure [TJ shows the ratio 
C[f(0)/\C[1\0)\ 2 versus T/d as a function of T/d on a 
semi-logarithmic plot for three cases as indicated in the 
figure caption. Case (i) with unequal transmission coef- 
ficients is obtained from experimental data and suggests 
that Eq. (fTB)) holds approximately for T > d. The two 



100 



* 'T—i 

u_ 



10 



Experimental data 



1 _ 



0.0 



0.4 0.8 

r/d 



1.2 



S io 

S3 



RMT simulation 



JeS . 



o 9_o.t>_e»e-S- e*> eP-oMi 



r/d 



Figure 1: Ratios of four-point- and squared two-point- functions of 
on a semi-logarithmic plot for three cases. Upper panel: Ex- 
perimental data 3, |l] (filled circles). Lower panel: Data generated 
numerically from Eq. JT} for A = 32 channels (open circles) and 
for A = 52 channels (full squares) with identical transmission coeffi- 
cients. 

The origin of the failure of Eq. (fTB)) for small values 
of T/d is easily understood. Qualitatively speaking, it is 
intuitively clear that cross-section fluctuations (measured 
in units of the average cross section) are much larger for 
isolated than for overlapping resonances. Quantitatively, 
the assumption that underlies Eq. (fTB)) is that the dis- 
tribution of the fluctuating S'-matrix elements S^ b (E) is 
Gaussian. That assumption holds only if the width 7 of 
the distribution is sufficiently small. Indeed, the unitarity 
condition Eq. (f6)) implies that for all a and b we must have 
|Sf fc | 2 < T a < 1, and 7 must be so small that the contri- 
bution of the tails of the distribution that extend beyond 
1/2 

the values ±T a , is negligible. Otherwise, significant de- 
viations from a Gaussian distribution are to be expected. 
To estimate 7 in the Ericson regime we use Eq. (fTT)) which 
for a b gives 7 = ^(O)! 1 / 2 = [T a T b / £ c TjVa. Thus, 
since ^2 C T C ^> 1 this estimate yields 7 <C 1 as expected. 
For smaller values of T/d, constraints on the distribution 
of S^ b due to the unitarity condition ([6]) are expected to 
become increasingly important as the number A of terms 
in the sum in Eq. (O decreases. 

We have checked this explanation by investigating the 
distribution of S^ 2 - A bivariate Gaussian distribution for 
Sf 2 implies that the distribution of z = \Sf 2 \/\Si 2 \ has the 



3 



form P(z) = (tv/2)z exp[— (7r/4)z 2 ], and that the phase 
of Sf 2 is distributed uniformly in the interval {0,27r}. In 
Fig. and Fig. [3] and for values of T/d indicated above 
each panel, we compare in the upper row the function P{z) 
with experimental and numerical data, respectively. In the 
panels in the lower row, we show the corresponding distri- 
butions of the phase. The numerical data in the three 
panels of Fig. [5] are obtained by simulating absorption in 
the resonator in terms of a large number of fictitious chan- 
nels with small transmission coefficients in each channel. 
Then their sum r is the only parameter. It was deter- 
mined as described in [3j from a fit of the experimental 
autocorrelation function to the analytic result given in |l| . 
The data show that agreement with the bivariate Gaus- 
sian distribution is attained for T/d > 1. The data in 
Fig. EH (generated numerically for 32 equal channels) show 
agreement with the bivariate Gaussian distribution only 
for T/d > 3. This suggests that the limit of a bivariate 
Gaussian distribution for Sz b with a ^ b is attained for 
larger values of T/d when all transmission coefficients are 
equal than when the transmission coefficients differ. These 
results account for the deviations from Eq. (fT3| displayed 
in Fig. □ 
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Figure 2: Distribution of the modulus (upper row) and phase (lower 
row) of Sf 2 obtained from the experimental data (histograms) and 
from numerical simulations (solid lines) for three values of T/d. In 
the upper row, the bivariate Gaussian distribution expected in the 
Ericson limit is shown as a dashed line. The transmission coefficients 
have the following values \3, .6]. For T/d = 0.17: Ti = 0.097, T 2 = 
0.130, r abs = 0.810, for T/d = 0.54: Tj = 0.417, T 2 = 0.475, r abs = 
2.274, for T/d = 1.01: T\ = 0.784, T 2 = 0.665, r abs = 4.903. 

Concerning the £-dependence of C}^ (e) and of | (e) | 2 , 
our data show that for those values of T/d where the ra- 
tio c[^(0)/\C^(0)\ 2 ~ 1, that dependence is sufficiently 
similar so that (e) can be reliably approximated by 
( e )| 2 - (By that we mean that the full two-point func- 
tion defined in Eq. (|4|) and not the approximate form 
Eq. (|T4"|) has to be used). That leaves us with the ques- 
tion how to approximate C[ 2 (s) analytically for T < d. 
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Figure 3: Distribution of the modulus (upper row) and phase (lower 
row) of 5^ obtained from numerical simulations (histograms) for 
A = 32 channels with identical transmission coefficients for three 
values of T/d. In the upper row, the bivariate Gaussian distribution 
expected in the Ericson limit is shown as a dashed line. 



An obvious possibility is offered by the function F^ 2 \e) 
defined in the first of Eqs. (fT2"]l . Figure [1] shows experi- 
mental (upper two panels) and numerically (lowest panel) 
generated values of \F± 2 and of c[^(s) versus e for 
two different values of T/d, and for A — 32 channels with 
identical transmission coefficients T a — 0.2 and a value of 
T/d close to the largest achieved in the experiment, re- 
spectively. The function \Fj£ (e) | was rescaled so that at 
e = it agrees with C[ 2 (0). In the upper two panels 
the experimental curves for c{ 2 ' (e) (filled circles) and the 
rescaled function | ^12^ ( e ) I are shown together with the an- 
alytic result for 1^12 ( e )l multiplied with the same scaling 
factor as the experimental one. In the lowest panel we 
compare the numerical result for (e) with the analytic 
rescaled result for \F^ 2 We note that the agreement is 
excellent for all three cases. Similarly good agreement was 
found also for A = 52 channels and several values of T/d. 
We conclude that we have reached a simple and reliable 
prescription for approximating C[ 2 (e) analytically for all 
values of T/d: Calculate cj^(0) and |-Fi2^( £ )| analytically. 
(The formulas needed 0, Q are given in the Appendix). 



Rescale \F^ 2 (e)\ so that its value at e = agrees with 
(7^2 (0). The resulting rescaled function l-F^O 5 )! is an ex- 
cellent approximation to c[ 2 \s) for all values of T/d. 

4. Elastic Case (a — b) 

The elastic case is analogous to the inelastic one only 
if S aa = 0. Then all conclusions drawn in Section [3] apply. 
In general, that is not the case. In particular, for T <C 
d (isolated resonances) the last term in Eq. ([T]) tends to 
zero and S aa approaches unity. The unitarity constraint 
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Figure 4: The dependence on e of C^(e) and of the renormalized 
function l-Fj^C 5 )! f° r three cases with r/d as indicated in the panels. 
In the upper two panels we show experimental curves for ell (s) 
(filled circles) and for l^^'C^)! (open squares), in the lowest panel a 
numerically generated curve for (e) (solid line) for 32 identical 

channels. In all three panels, the analytic function |i*i2^(£)| i s shown 
as dashed line. Both the experimental and the analytic values for 
| '(e) I are renormalized with a factor determined from the analytic 
value for c[f(0). 

on S implies \S aa + S^ a \ < 1 and the distribution of S^ a 
must then become skewed. Therefore, we have to expect 
that Eq. ([T3"]l applies less generally in the elastic than in 
the inelastic case, and that C^(e) in Eq. Q plays an 
important role. Actually, it was pointed out in Refs. [J, 
[B[ that even in the Ericson limit, a profound difference 
between the elastic and the inelastic cases exists: For T ^> 
d, \Sa a \ 3 and \S^ a \ 4 have similar values unless T a w 1 or 
Saa ~ 0. That shows that in order to predict C aa {e) we 
need to know both d 3 J (e) and Caa (e) for ah values of 
T/d. 

Concerning c[i (0), we proceed as in Section [3] and 

display in Fig. [5] the ratio C[f(0)/\c[f(0)\ 2 in the upper 
panel for experimental data, in the lower one for a numeri- 
cal simulation with A = 32 equal transmission coefficients, 
as indicated in the figure caption. We reach the same con- 
clusions as for a ^ b: When the transmission coefficients 
differ, the ratio is close to unity for r > d, when they are 



equal for T > 3d. As in the inelastic case, the reason for 
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Figure 5: Same as Fig. ^ but for S 1 ^. 



the failure of Eq. (|13[) is the non-Gaussian distribution of 
S & . This is shown in Figs. [6j [71 The deviations from a 
bivariate Gaussian distribution now occur for larger val- 
ues of T/d than in the inelastic case, in keeping with the 
expectations formulated at the beginning of this Section. 
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Figure 6: Same as Fig. [2] but for S^. 

In Fig. [S]we show that, as in the elastic case, the rescaled 
function |-Fn^(£)| agrees well with cj^(e) for all values of 
T/d. Such agreement was likewise found for the case of 
A = 52 channels with identical transmission coefficients. 
Figure [5] shows that the three-point function 9le \ C[^(s)J 

is approximated quite well by the (unrenormalized) func- 
tion SKe [Fix (e)j defined in the second of Eqs. p2[l for the 
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Figure 7: Same as Fig. [3] but for 5^. 

experimental data with T/d = 0.17 (upper panel), where 
the resonances are nearly isolated, and for the numerical 
simulations with A = 32 identical transmission coefficients 
T a = 0.2 (lower panel). However, considerable deviations 
are observed for the experimental data with T/d = 1.01 
(middle panel). Thus, the conclusions drawn for the inelas- 
tic case apply similarly to the elastic one only in the regime 
of isolated resonances and the Ericson regime, where S aa 
is vanishingly small. To obtain an analytic expression for 
the cross-section autocorrelation function we now need to 
replace G$ (e) with the rescaled function |F-^ (e) | and the 

function me (c[f (sfj with «e (F^(e)\ While the for- 
mer replacement yields an excellent approximation for all 
values of T/d, this is not generally true for the latter in 
the intermediate regime of weakly overlapping resonances. 
Still, for the numerical simulations with A = 32 equal 
transmission coefficients T a = 0.2, which corresponds to 
T/d = 1.02, the resulting analytic expression for the cross- 
section autocorrelation function yields a good approxima- 
tion, because there both, the values of IHe (F^iej) and of 

d\e (c^(e)\ are negligibly small as compared to C a &(e), 
whereas this is not true for the experimental data with 
T/d =1.01. 

5. Summary 

We return to the questions raised in the Introduction. 
We have shown that Eq. (fT3"]) holds essentially only in the 
Ericson regime and, for the elastic case, even there only 
approximately. This is because of the constraint imposed 
on the distribution of S R by unitarity. Hence, in the inelas- 
tic case the cross-section autocorrelation function can be 
reliably predicted from Eq. (THJ|) for T/d > 1 for incquiva- 
lent channels and for T/d > 3 for identical channels. Even 
more stringent constraints exist in the elastic case. We 
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Figure 8: Same as Fig. [4] but for the elastic case. 

have demonstrated, however, that for the inelastic case 
for all values of T/d and for the elastic case under certain 
conditions an excellent approximation to the cross-section 
autocorrelation function may be obtained by replacing in 
Eq. © the three-point function C$(e) by F^\s) and 
the four-point function Cj$(e) by the rescaled function 
F^t\e). With these replacements, Eq. ^ takes the form 



C ab (e) 



25 ab Vttl S, 



Saa F<$(e) + S aa Fg\-e) 



^(0)1 



(15) 



The input parameters for the evaluation of Eq. ([15)) are the 
average i'-matrix elements and the associated transmis- 
sion coefficients for all channels and, for the dependence on 
e, the average level spacing d of the scattering system. In 
terms of these parameters, the two-point function cf? (e) 
is given in Ref. [j| , and (0) and the functions (e) 
and F^\e) in Refs. 0, For the convenience of the 
reader, all relevant formulas are collected in the Appendix 
in a form suitable for numerical implementation. For the 
elastic case, Eq. (1151) provides an excellent approximation 
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Figure 9: The dependence on e of 9te fi^f (e)) and of <Ke (c{f (e) 
for three cases with F/d as indicated in the panels. In the upper two 
panels we show experimental curves for <Ke (c[f(ef\ (filled circles) 

and for ?Re M*i?( e )j (open squares), in the lowest panel a numer- 
ically generated curve 9te (c^f 1 (e) J for 32 identical channels (solid 

line). The analytic function 5Ke (F^(e) \ is shown as dashed line in 
all three panels. The limited accuracy of our comparison is indicated 
by the difference between the dashed line and the open squares in 
the upper two panels. 



for the cross-section autocorrelation function for all values 
of T/d only for e = 0. For e > Eq. (fT5]) yields a good 
approximation in the regime of isolated and strongly over- 
lapping resonances, whereas in the intermediate regime of 
weakly overlapping resonances this is true, if the contri- 

(3) (3) 

bution of C ab (e) or equivalently of \F^ b '(e)\ to C a b{s) is 
negligible. 
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Appendix 

For the evaluation of C a t>( £ — 0) in terms of S R and S as 
given in Eq. © we need to compute C^ b (e = 0), C a 3 J(s = 



0) = ^(0)|S« fc (0)| 2 andCf ( £ = 0) = |S a fl b (0)lH^(0)i 2 
The autocorrelation coefficient C^ b \e = 0), S^(0) \S% b (0)\ 2 
and |52»,[0)| 4 are given in terms of a threefold integral, c.f. 
Refs. [U 



^(0)= (16) 

-i pOC pOC pi 

- / d\i / dA 2 / d\J(X,X 1 ,X 2 ) 
8 Jo Jo Jo 



1-T,A 



r%\x,xi,x 2 ) 



V(l + T c Ai)(l+T c A 2 ) 
where the integration measure is given as 



J(X,Xx,X 2 ) 



(17) 



A(1-A)|A!-A 2 | 



(A + Ai) 2 (A + A 2 ) 2 v / Ai(l + Ai)A 2 (l + A 2 ) 



For n = 2, i.e. for |Sf 6 (0) | 2 the factor JF(A, Ai, A 2 ) equals 



"^('V Ai, A 2 ) - 

Sab \S aa \ T a 



2 rl ^1 ^2 2A 



(1 + 5 ab ) T a T b 



1 + T a X 1 1 + T Q A 2 1 - T a X 
Ai(l + Ai) 



(18) 

2 



(l + TaXJil + nXt) 

A 2 (l + A 2 ) 2A(1-A) 



(l + T a A 2 )(l + T fc A 2 ) (l-T a A)(l-T b A) 



for n = 3, i.e. for (0)|5f b (0)| 2 it equals 

^ ) (A,A 1 ,A 2 )= (19) 
- 5"aa ^4trg(/i a i/ a ) -I- 2trg(/x a )trg(z/ a ) 

+ r„|trg (p 2 ) +i[trg(^ Q )] 2 }|^ Qb 



and for rt = 4, i.e. for 1 S^ b (0)\ 4 we have 
•^ib ) (^' ^i; ^2) = 

1 2 

trg(y a ) + -r a [trg(/i )j 



/ <j trg(i^) + ^r fc [trg(^ 6 )] 2 



(20) 



+ 6 ab ( [trgK)] 2 + 4trg (z/ 2 ) 

+ r 2 trg> 2 ) {trg (^ 2 ) + [trg(/x a )] 2 } 
+ r a | [trg(|U a )] 2 trg(i/ tt ) 

+ 8 trg(Mo)trg(/Xo«/o) + 8 trg {[i\v a ) } 



Here, r a = 1 — T a and fi a , v a are the matrices 

[i a = T a X (1 + TqAo) -1 , 

v a = T a 2 A (l + A )(l+T a Ao)" 2 , 



(21) 



where Ao is the 4x4 diagonal matrix with entries Ax, A2, 
—A and —A. For a 4 x 4 matrix M with diagonal ele- 
ments mu the graded trace trg is defined as trg(M) = 
(mix + 77122) — ( TO 33 + m 4i)- Analytic expressions for 
the functions F^\e) and F^(e) defined in Eqs. (fT2|) are 

( 2) 

obtained by multiplying J 7 ^ (A, Ai, A 2 ) in Eq. (fl"8j) and 

^i ) (A,Ai,A 2 ) in Eq. (59]) and .7^ (A, Ax, A2) in Eq. ® 
with the exponential 



exp [-i (Ax + A 2 + 2A) 7re/d] 



(22) 



The integral Eq. (|16p contains several singularities. These 
can be handled by proceeding as in Section 5 of Ref. 14 1 . 
We define the variable 



p = Ax + A 2 + 2A . (23) 

Then the exponential in Eq. (f2"2"|) turns into 

exp(— ipire/d) . (24) 

We distinguish the two cases p < 2 and p > 2. For p < 2 
we define two further integration variables 



a/Ax + A 



t = 



a/Ax + A 2 + 2A 

VAT 



VAx + A 

For the inverse transformation that yields 
A = ps 2 (1 - t 2 ) , 



Ax = ps 2 t 2 , 



A, 



2,2 



/j(l + .s z t 



2s 2 



(25) 
(26) 



(27) 
(28) 
(29) 



The Jacobian for this transformation is equal to 4p 2 s 3 t. 
The threefold integral in Eq. (TIB]) becomes 



F ab\ £ )(p<2) - 

r 2 rl/V2 rl 

dp I ds I dt 



(30) 



-ipire I d 



Jo Jo ~p(l-a 2 )V(l+Ps 2 * 2 ) 

[1-ps 2 (1-i 2 )] (l-2s 2 ) (1-t 2 ) 
^(1 + s 2 < 2 - 2s 2 ) [1 + p (1 + s 2 < 2 -^ 2 )] 



n 



1 -r r A 



^ ) (A,A 1 ,A 2 ). 



V /(1 + T C A 1 )(1+T C A 2 ) 
For p > 2 we define additional integration variables 

(31) 



s = vaT+a . 

AT 



t 



\Ai + A 



(32) 



For the inverse transformation that yields 
A = 



s 2 (l-< 2 ) 



Ax = s¥, 



A, 



p + s 2 t 2 - 2s 2 . 



(33) 
(34) 
(35) 



The Jacobian equals 4s 3 i and the threefold integral in 
Eq. (TTd) becomes 



F ab\ £ )(p>2) 



(36) 



00 ry/p/ 2 r 1 e —ipTre/d 
dp ds dt 2 

2 JO J Vl-s- 2 8(s-l) (p — S 2 ) 



X 



[l-s 2 (l-£ 2 )] (p-2s 2 ) (1-t 2 ) 
y/(p + sH 2 - 2s 2 ) (1 + sH 2 ) (1 + p + sH 2 - 2s 2 ) 

1 - T c A ^,(„) 



V(l + T c Ax)(l+T c A 2 ) 
where the (9-function is defined by 

9(x) = for x < , 
6»(.t) = 1 for x > . 
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